Topological superfluids with time reversal symmetry 
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It is shown that superfluids in two and three dimensions which have time reversal invariant ground 
states have phases which are distinguished by a topological invariant. Further, it is shown that the 
B-phase of 3 He is a superfluid in the non-trivial topological class. Superfluids in the non-trivial 
topological class are shown to have gapless edge states and support various kinds of vortices with 
zero energy modes localized in their cores. Some of these vortices have non-abelian statistics. 

PACS numbers: 



There has been a considerable amount of interest in 
phases characterized by a topological invariant (topo- 
logical phases). The most well known examples of such 
phases are the quantum Hall states 0,0, 03]- Chiral super- 
conductors and superfluids are another class of systems 
which are characterized by a topological invariant Q • A 
droplet of 3 He in the A-phase [5| is an example of a 2D 
chiral superfluid while strontium ruthenate is thought 
to be a chiral superconductor jf|. Chiral superconduc- 
tors and superfluids have a number of interesting prop- 
erties. A neutral chiral superfluid has a macroscopic cur- 
rent which runs along the edge of a bounded sample. 
In a superconductor, these currents are screened but still 
substantial @,S]- Chiral superconductors and superfluids 
have gapless edge states. They also support exotic defects 
called half quantum vortices which carry half a quantum 
of flux. Each half quantum vortex carries a Majorana 
fermion and a pair of Majorana fcrmions is equivalent 
to a Dirac fermion. A system of 2n vortices thus has a 
2 n dimensional vector space of degenerate ground states. 
The zero energy Majorana fermions are stable since they 
cannot couple to local operators. When the vortices are 
braided around each other, the vector corresponding to 
the ground state is rotated in the 2™ dimensional vec- 
tor space and t his p henomenon is known as non-abelian 
statistics [l,[l(l[ll|. 

Chiral superconductors and the quantum Hall states 
break time reversal symmetry. Until recently, examples 
of known topological phases were restricted mostly to sys- 
tems whose ground states broke time reversal symmetry. 
In this paper we study superfluids with unbroken trans- 
lational symmetry whose ground states are time reversal 
invariant in two and three dimensions. We find that the 
ground states of these systems in two and three dimen- 
sions can be classified by means of a topological invariant 
and that superfluids in the non-trivial topological class 
have a number of interesting properties. These systems 
have robust gapless edge states which are stable against 
perturbations that do not break time reversal symmetry. 
They can also support exotic defects which have zero- 
energy Majorana fermions and can lead to non-abelian 
statistics. 

There has been considerable recent work on materials 



and systems with time reversal symmetry which have a 
non-trivial topological Zi invariant and some of these are 
believed to have been experimentally detected [HI, [l3j ]. 
These systems are generically called quantum spin Hall 
systems [13]. We show that the B phase of superfluid 
Helium-3 can be identified as a "topological superfluid" 
in the non-trivial topological class. 

For the rest of this work, unless explicitly stated, by 
superfluids, we mean superfluids with unbroken transla- 
tion symmetry. Lattice superconductors which are not 
invariant under the full translational group have previ- 
ously been studied in an earlier work [15|. In Sec. I, we 
study and obtain a topological classification of superflu- 
ids in two and three dimensions. In Sec. II, we study the 
edge states in these superfluids. In Sec. Ill, we present ex- 
amples of superfluids in the non-trivial topological class. 
In Sec. IV, we consider exotic vortices and their statistics 
in these systems. 



I. TOPOLOGICAL SUPERFLUIDS WITH TIME 
REVERSAL SYMMETRY 

The mean field BdG Hamiltonian which characterizes 
a superfluid may be written in the form: 



H 



where 



BdG 



H 
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h(k) A(fe) 
At(Jfe) -h T (-k) 



(1) 
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Here h represents the single particle Hamiltonian, A 
the order parameter that characterizes superconductiv- 
ity and i/j,ip are the fermionic operators in the two 
component Nambu formalism. The BCS ground state 
wavefunction of the superfluid is annihilated by opera- 
tors of the form j(k) = J2 a =t,l u k, a tpk, a + v k , a T_ kja 
where (uk, a , Vk, a ) T is a negative energy eigenvector of 
the matrix H. 

As long as the system has a bulk gap, at each point 
in momentum space, there are two eigenvectors of this 
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matrix, say ei(fe) and e2(fe) with negative eigenvalues. 
As k — > oo, the two dimensional vector space spanned 
by these states, which we denote by V(k) goes to a fixed 
two dimensional space 2l| . As far as the topology of the 
ground state wavefunction is concerned, the base space 
which is momentum space is thus a sphere, S n obtained 
by the one point compactification of R n , where n is equal 
to 2 or 3 in the cases of our interest. 

When the ground state wavefunctions has time reversal 
symmetry, the eigenvectors of H at k,~k are not inde- 
pendent. Time reversal symmetry requires that if u is an 
eigenstate at fc, then Qu is also an eigenstate with the 
same energy at —k. While the two dimensional vector 
space consisting of the negative energy eigenvectors can 
be written in terms of a continuous basis locally in mo- 
mentum space, finding a global continuous basis is not 
always possible. The ground state thus defines a twisted 
two dimensional vector bundle on S 2 or S 3 as the case 
may be. 



Superfluids in 2D 

Let us first study the topology of the ground state 
wavefunction of a two dimensional superfluid follow- 
ing We divide S 2 into two patches, E and E' where 
E = {k 3 k < 2} and E' = {k 3 k > 0.5}. On each of 
these patches we define a basis of two continuous orthog- 
onal vector functions, |ui,»(Ab)), |«2,i(k)) & V(k) where 
the index i is a label for the patch E or E'. We repre- 
sent the vector Uj(fe) = c\Vi^(k) + c 2 v 2 ^(k) as the spinor 
(ci,C2) T . The operation of time reversal in this basis 
then corresponds to the operator ia^K^ where Kq is 
the complex conjugation operator and <I> is the involu- 
tion operator that takes vectors at k to vectors at — k. If 
a vector v is represented by the spinors ve and ve> in the 
two different bases defined above, then ve and ve> can be 
related by means of a transformation matrix U such that 
ve = Uve> where U is an element of the unitary group, 
U(2). The transition function, U(k), which glues the vec- 
tor representations VE(k),VE'(k) on the two patches, on 
the circle T — {k 3 k = 1} € Ef] E' can be written as 
e ia{k)i e in(<j>).<Tu{k) _ ^hen time reversal invariance leads 
to the condition: 



U(-k) 



-ia(k)I ^in((p) .eruj(k) 



(3) 



Then, as shown in previous work 151 ] . U2D, the set 
of all transition functions for these ground states in two 
dimensions consists of two topologically distinct classes 
of transition functions. The trivial class corresponds to 
the case when it is possible to find a globally continuous 
basis for the vector space spanned by ei(fc), e2{k). The 
non trivial Z 2 phase corresponds to the case when the 
transition function cannot continuously be deformed to 
unity. A simple formula for the Z2 invariant can be writ- 
ten down when the wavefunctions ei(fc),e2(fc) have well 



defined Chern numbers. When this is not the case, one 
can still always find a set of two orthogonal "wavefunc- 
tions" e[(k),e 2 (k) 6 V(k) for which well defined Chern 
numbers exist. The Z2 invariant can then be expressed in 
a particularly simple form in terms of these Chern num- 
bers. It is given by |c„| mod 2 where c n is the Chern 
number of either of these functions. 



Superfluids in 3D 

Now consider a superfluid in three dimensions. We 
again divide momentum space, which is now topologi- 
cally equivalent to the three sphere S 3 for our purposes, 
into two patches E = {k 3 k < 2} and E' = {k 3 
k > 0.5}. As before, on each of these patches we define 
a basis of two continuous orthogonal vector functions, 
|«i,»(fe)), \v2,i{k)) € V(k) where the index i is a label for 
the patch E or E\ 

The transition function, U(k) which glues the vec- 
tors on the two patches, VE(k),VE'{k) on the two sphere 
T = {k 3 k = 1} e Ef]E' is again an element of U(2) 
and can be written as e la ( fc ) / e l ™W)-<M fc ). Time rever- 
sal symmetry results again in the constraint given by 
Eq. ©■ Let C be any great circle on T. Then the func- 
tion obtained by restricting the domain of U (k) to T is 
a function which belongs to the class of two dimensional 
transition functions U2D- This thus defines a projection 
from the class of transition functions on S 2 for the three 
dimensional case, which we call U3D, to the class of tran- 
sition functions for the two dimensional case, U2D- As 
previously discussed, there are two topologically inequiv- 
alent classes in U2D- 

Further, it is well known that tt 2 ([/(2)) = tt 2 (SU(2)) 
I ), This implies that the class of transition functions in 3D 
is determined by the topological class of the 2D transition 
function that it projects to. 

Thus the topological invariant for the 3D superfluid 
may also be evaluated using an extension of the simple 
Chern number formula from the 2D case. Let P be any 
two dimensional plane which maps onto itself under time 
reversal and contains the origin and let e' 1 (k),e 2 (k) be 
a set of orthogonal wavefunctions which span V(k) at 
each point k in the plane, which map onto each other 
under time reversal symmetry and for which well defined 
Chern numbers exist. Then, the 3D Z2 invariant is sim- 
ply |c„| mod 2, where c n is the Chern number of one of 
the wavefunctions. 

We have shown above that there are thus two distinct 
topological phases of superfluids at the mean field level. 
This implies that if the full Hamiltonian of the system 
is adiabatically perturbed, as long as the ground state 
does not spontaneously break time reversal symmetry, 
the ground state remains in one of two distinct phases, 
unless there is a phase transition. 

We note that superfluids in 2D have a single Z 2 invari- 
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ant as do lattice superconductors. However in 3D, lat- 
tice superconductors have four Z2 invariants 15], while 
superfluids with unbroken translational symmetry have 
a single Z 2 invariant. This invariant corresponds to the 
fourth Z2 invariant which is intrinsically three dimen- 
sional in nature and corresponds to the invariant in in- 
sulators and lattice superconductors which determines 
whether the system is in the strong or weak topological 
class. 



II. EDGE STATES IN TOPOLOGICAL 
SUPERFLUIDS 

Superfluids in 2D 



Hamiltonian H(ak x + (3k y = 0, z) represents a 2D su- 
perfluid Hamiltonian which from our previous discussion 
is a 2D TR invariant Hamiltonian in the non-trivial Z 2 
class. It therefore follows that the eigenstates of this 
Hamiltonian has a set of robust edge states with a ID 
Dirac spectrum. Since this is true for arbitrary a, f3, it 
follows that the eigenstates of the Hamiltonian in general 
have a 2D Dirac energy spectrum at low energies. While 
the spectrum is Dirac-like, the eigenstates are Majorana 
rather than Dirac fermions. There is a single pair of zero 
energy Majorana fermions at the edge corresponding to 
the k x — k y — state. More generally when there are 
open boundary conditions, there are two surfaces. From 
the above argument, it is expected that each one will 
have a single pair of zero energy Majorana fermions [22] ]. 



A chiral 2D superconductor and a quantum Hall in- 
sulator are both characterized by the same topological 
number, namely the first Chern number The con- 

nection between the Chern number and the presence of 
robust edge states is well known in the context of the inte- 
ger quantum Hall effect 17, 18j]. Further, it is also known 
that these edge states exist in chiral superconductors and 
superfluids. This follows both from the Chern number - 
edge state connection and also from explicit calculations 
using the BdG Hamiltonian 0]. The gapless edge exci- 
tations in the case of a chiral superfluid/superconductor 
are Majorana fermions or equal linear combinations of 
particles and holes. In the case of insulators with time 
reversal symmetry, the Z2 invariant determines whether 
robust edge states exist or not. A superconductor whose 
ground state respects time reversal symmetry can not 
have net charge currents along the edge. Nevertheless 
since the ground state of a non-trivial 2D Zi supercon- 
ductor may be adiabatically continued to a product state 
of two wavefunctions each of which has the same form in 
position space as the ground state wavefunction of a chi- 
ral superconductor with an odd Chern number, it follows 
that the TR invariant superconductor of the non trivial 
Z2 class has a robust pair of edge states. When the Chern 
number is one, there is precisely one set of edge states 
and a single pair of Majorana fermions at zero energy 
at the edge. This pair of Majorana edge states is stable 
as long as time reversal symmetry is preserved because 
no operator that is invariant under time reversal symme- 
try can have a non-zero matrix element between the two 
states fiHliol]. 



Superfluids in 3D 

We now study the 3D case. Consider a semi infinite 
three dimensional superfluid with periodic boundary con- 
ditions in any two directions, say in the x and y directions 
and which exists in the region z < 0. The momentum 
variables k x ,k y are then good quantum numbers. The 



III. EXAMPLES OF TOPOLOGICAL 
SUPERFLUIDS 



In this section, we shall present some examples of states 
in two and three dimensions that are in the non-trivial 
Z2 class. Since these states are triplet states, it is useful 
to introduce the d-vector notation. The order parame- 
ter of a general triplet superconductor/superfluid can be 
written in the form [5j: 



A = A Q I id », dx ., dz , 
d z idy + d x 



(4) 



where in the absence of textures, defects and edges, d is 
a function only of k. 

The state whose d-vector is given by 

d(k) = i(k x + ik v )y 

is an example of a chiral p+ip superconductor/superfluid. 
On the other hand, the state with 



d(k) 



Ky 



(5) 



is an example of a superfluid with a non-trivial Z2 in- 
variant [23] ]. Since the spins decouple for such a state, 
the up and the down spins may be considered separately. 
The wavefunction for the up/down-spins are identical to 
those of a chiral p x — ip y /p x + ip y superfluid. Thus, the 
form of the wavefunctions and their spectrum for chiral 
superfluids in rectangular and circular geometries with 
an edge may be used to deduce the wavefunction for 
the ground state of the above system, but shall not be 
presented explicitly here. It follows from the analogy 
with chiral superfluids, that there is a single pair of sta- 
ble zero energy Majorana edge states as previously stated 
in Sec. II. 

The Balian Werthamer state which is believed to exist 
in the B-phase of 3 He has the description : 



d{k) = k 



(6) 
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This is a gapped 3D superfluid which is time reversal in- 
variant. To determine its topological invariant, we con- 
sider the plane k z = in momentum space. The order 
parameter then corresponds to the 2D superfluid given 
by Eq. © which is in the non-trivial Z2 class of 2D su- 
pcfluids. From the results of Sec. I, it follows therefore 
that 3 He - B is a superfluid in the non-trivial Z2 class. 

Indeed, as anticipated in Sec. II, the superfluid state 
when confined to a boundary has gapless edge states. 
Since the BdG Hamiltonian for this state is invariant un- 
der simultaneous rotations in spin and orbital space, one 
may consider without loss of generality a three dimen- 
sional sample confined with an infinite wall at y=0. We 
further consider the set of states with k z = 0. The re- 
sulting system can be regarded as the sum of two two 
dimensional chiral superfluids with opposite spin and chi- 
rality (24|. The edge states of a chiral 2D superfluid in 
a rectangular geometry were analyzed in Ref . [7fl . It was 
found that the spectrum was linear in the momentum 
component parallel to the edge and that there was pre- 
cisely one zero energy Majorana mode. It follows that 
the confined 3D superfluid considered here has a spec- 
trum E oc yfel + fcf which is 2D Dirac-like. 



These are spin- vortices which may be seen as arising from 
the slow rotation of the d(k) in position space. The over- 
all phase of the order parameter does not vary and thus 
these defects preserve time reversal symmetry. The spins 
are still decoupled and each spin component can be re- 
garded as a chiral superfluid which now has a full quan- 
tum vortex but of opposite relative chirality. Further, 
it follows that there is a single pair of zero energy Ma- 
jorana fermions which are stable to local perturbations 
which do not break time reversal symmetry. The statis- 
tics of these fermions for each individual spin component 
is non-abelian 25| and this could possibly be harnessed 
for quantum computation if the spins are kept decoupled. 

The 3D superfluid previously discussed, which is de- 
scribed by the equation: d(k) = k is also degenerate 
with the set of states d(k) = R[k] where R is any ro- 
tation in real space when the system has both spin and 
orbital rotational symmetry. Since R is an element of 
50(3) and 7Ti(50)(3) = Z 2 , the superfluid then sup- 
ports line defects which correspond to non-contractible 
paths in 50(3). An example of such a vortex is 



d{k,0) = R(n,e)[k] 



(10) 



IV. EXOTIC VORTICES AND NON-ABELIAN 
STATISTICS 



The 2D topological superfluid state whose order pa- 
rameter is given by Eq. (0 supports many kinds of exotic 
vortices. In the vortices that we shall be interested in, 
the order parameter takes the form: 



A(r,0) = A (r) 



(I,. 







Here r, 9 are the polar coordinates in the vortex and <f>± 
are the phase windings of the different components of the 
order parameter. When <f> + = 9,4>- = or vice versa, 
the defect is a half quantum vortex. In a superconduc- 
tor such a defect would carry half a quantum of flux. 
Only one of the spin components are involved in the low 
energy physics. We may therefore consider a spinless chi- 
ral superconductor with a full quantum vortex, which as 
pointed out in Ref. Q has a single zero energy Majorana 
core state at the center of a vortex. Further the statistics 
of these states is non-abelian. 

When the system has both spin and orbital rotational 
symmetry the state given by Eq. ([5]) is degenerate with 
the set of states 

d(k) = R[k] (8) 

where R[k] represents the vector obtained by a rotation, 
R, of the k-vector in the two dimensional momentum 
space plane. One can therefore consider defects which 
are represented by 



±0 , 4>_ = =p0 



(9) 



where 9 is the polar coordinate in the vortex and 
R(n,9)[k] is the vector obtained by rotating k through 
the angle 9 about the axis n. To determine the low en- 
ergy states, we apply periodic boundary conditions in 
the z direction and set n = z. Then for the set of states 
k z = 0, we see that the order parameter at the vortex 
has the same form as the defect in the 2D superfluid con- 
sidered in Eqs. and (|9|). There is thus a stable pair of 
zero energy Majorana states at the core of these vortices. 
The 3D superfluid state also hosts analogs of the half 
quantum vortex. The discussion from the 2D case can 
be easily carried over to the defects in the 3D superfluid 
states. 

The above discussion pertained to superfluids with 
unbroken translational symmetry. Lattice supercon- 
ductors in two dimensions are described by a single Z2 
invariant, while in 3D, they are described by four such 
invariants. The fourth Z2 invariant for superconductors 
is analogous to the single Z2 invariant for superfluids 
discussed above. The discussions pertaining to the edge 
states and the exotic defects presented above is also valid 
for superconductors with a non-trivial Z2 invariant in 2D 
and a non trivial fourth Z2 invariant in three dimensions. 

In conclusion, we have shown that superfluids which do 
not break translational symmetry have two phases deter- 
mined by topology. We also showed that the non-trivial 
topological class of superfluids have a number of inter- 
esting properties such as edge states and exotic vortices 
which obey non-abelian statistics. We also provided ex- 
amples of such superfluids and showed that the B-phase 
of 3 He is an example of a non-trivial topological super- 
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fluid. It is hoped that the present work will spur further 
interest and work in these systems. 
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